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ABSTRACT 



< 

i-C ' We present results from a wide-area photometric survey of the Phoenix dwarf 

galaxy, one of the rare dwarf irregular/ dwarf spheroidal transition type galaxies (dTs) 
q ' of the Local Group (LG). These objects offer the opportunity to study the existence 

?-h , of possible evolutionary links between the late- and early- type LG dwarf galaxies, 

c/3 ■ since the properties of dTs suggest that they may be dwarf irregulars in the process 

of transforming into dwarf spheroidals. 

Using FORS at the VLT we have acquired VI photometry of Phoenix. The data 
reach a S/N~10 just below the horizontal branch of the system and consist of a mosaic 
| of images that covers an area of 26' x 26'centered on the coordinates of the optical 

■^j- • center of the galaxy. 

' Examination of the colour-magnitude diagram and luminosity function revealed 

t^-* , the presence of a bump above the red clump, consistent with being a red giant branch 

■ bump. 

' The deep photometry combined with the large area covered allows us to put on 

, a secure ground the determination of the overall structural properties of the galaxy 

£SJ ' and to derive the spatial distribution of stars in different evolutionary phases and age 

ranges, from 0.1 Gyr to the oldest stars. The best-fitting profile to the overall stellar 
population is a Sersic profile of Sersic radius Rg =1.82'±0.06'and m =0.83±0.03. 
We confirm that the spatial distribution of stars is found to become more and 
>^ \ more centrally concentrated the younger the stellar population, as reported in previous 

5— i ■ studies. This is similar to the stellar population gradients found for close-by Milky Way 

dwarf spheroidal galaxies. We quantify such spatial variations by analyzing the surface 
number density profiles of stellar populations in different age ranges; the parameters of 
the best-fitting profiles are derived, and these can provide useful constraints to models 
exploring the evolution of dwarf galaxies in terms of their star formation. 

The disk-like distribution previously found in the central regions in Phoenix ap- 
pears to be present mainly among stars younger than 1 Gyr, and absent for the stars 
> 5 Gyr old, which on the other hand show a regular distribution also in the center of 
the galaxy. This argues against a disk-halo structure of the type found in large spirals 
such as the Milky Way. 

Key words: techniques: photometric - Local Group - galaxies: individual: Phoenix 
- galaxies: stellar content - galaxies: structure - galaxies: evolution 



* Based on FORS observations collected at the ESO, proposal 
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1 INTRODUCTION 

Dwarf galaxies are the most common type of galax ies in the 
nearby universe (see e.g. iKarachentsev et al.|[2004 ). and are 
expected to be among the first systems to form in the cos- 
mological context. At larger redshifts galaxies such as those 
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typically seen in the Local Group (LG) cannot be observed. 
Therefore detailed observations of current properties of LG 
galaxies offer the opportunity to gain insight into the earli- 
est stages of galaxy formation in the most numerous type of 
galaxies in the Universe. 

LG dwarf galaxies are classified in two main categories: 
the late-type objects, containing gas and currently form- 
ing stars, often have irregular appearance in the optical and 
therefore are called dwarf irregulars (dlrr), while the early- 
type dwarf spheroidals (dSph) - devoid of neutral gas - show 
a m ore regular m orphology and no current star formation 
(e.g. lMateol["l998l ). Dwarf galaxies with intermediate prop- 
erties - such as no current star formation but containing gas 
- are called transition types (dTs). 

An open question is whether the various types of dwarf 
galaxies in the LG are intrinsically different objects, or 
whether they descend from the same progenitors and have 
evolved through different paths because of environmental 
and/or internal processes. 

The continuum of p roperties shown by t he various types 
of LG dwarf galaxies ijTolstov et all 12009"). which makes 
their classification sometimes uncertain, seems to argue in 
favour of an evolutionary link. For example, PegDIG is 
the most luminous of the few LG transition type dwarfs 
(Phoenix, LGS3, DDO210, Leo T, VV 124) and, because 
of this, it is sometimes considered as a dlrr. For Fornax, 
which has formed stars un til very recently (50-100 Myr ago 
Coleman fc de Jondl200"&T ) but is now devoid of gas, its clas- 
sification as a dSph seems to be merely the result of the 
particular moment in time in which we are observing this 
galaxy: if we had been able to observe Fornax a few hundreds 
Myr ago, while still forming stars, we would have classified 
it as a dT or even as a dlrr. 

Another hint to a possible evolutionary link is the ex- 
istence of a "morphology-density" relation: dlrrs are found 
at relatively large distances from the large LG spirals (i.e. 
> 300 kpc from the Milky Way and M31), while the great 
majority of early type dwarf galaxies are satellites of the 
Milky Way or M31, being located at < 300 kpc from them. 
This already suggests that the interaction with the large 
spirals (hereafter referred to as "environment") must have 
played a role in determining the different evolutionary paths 
of late and early type dwarfs. Models based on N-body sim- 
ulations do show that such interactions are a viable expla- 
nation to the obse r ved "morphology-densit y" relation (e.g. 
iMaver et ail 120061 : iKazantzidis et al.l [20TH ) . However, the 
presence of dSphs (e.g. Cetus and Tucana) and dTs (e.g. 
Aquarius) at distances such that any strong interaction with 
either the MW or M31 can be excluded indicates that the 
environment cannot be the only factor at play but that also 
internal factors (e.g. the depth of the potential well of the 
single objects) might play a role. 

The lack of knowledge about the detailed properties of 
dlrrs and dTs hinders the exploration of a possible common 
origin of LG dwarfs, and what might have caused a different 
evolutionary path. 

At a distance of 415 kpc l|Hidalgo et al. 2009, hereafter 
H09) , Phoenix is the closest transition type dwarf, but even 
basic pro perties such as its e x tent a re rather uncertain. For 
example, Ivan de Rvdt et al.l (Il99ll, hereafter VDK91) de- 
rived an extent of 8.7', while iMartmez-Delgado et alj (| 19991 . 
hereafter MD99) found a larger value, with a nominal tidal 



radius of r t = 15.8_ 2 g ar cmin. Also its heliocentr ic systemic 
velocity is under debate: llrwin fc Tolstoy (12002 ) der ived a 
velocity of -13 ±9 km s" 1 . while IGallart et all feoOll ) mea- 
sured a velocity of — 52±6 km s _1 . Both studies though con- 
sider likely an association between Phoenix and the cloud 
of H I gas found around the g alaxy at heliocentr ic velocity 
between -30 and -14 km s" 1 l|Young et al.ll2007l '). 

The star formation history (SFH ) of Phoenix has 
inste ad been exte n sively studied (see van de Rvdt et ail 
I l99ll: lHeidet al.1 Il999l; iMartmez-pelgado et all Il999l; 



i Holtzman et all l2000l: IGallart et al.l 120041 : iMenzies et all 
20081: iHidalgo et all l2009h . The co-existence of classical 
Cepheids and anomalous Cepheids and RR Lyrae variables 
in this system a rgues for an extende d and complex star for- 
mation history (|Gallart e t al. 2004). Phoenix has a similar 
luminosity to the Sculptor dSph (Ly = 0.9 x 1O 6 L0 and 
L v = 2.15 x 10 6 L Q , respectively; Mateo 1998), and like 
Sculptor it formed most of its stars more than 10 Gyr ago 
(H09). However, unlike Sculptor or the majority of other 
dSphs it does show recent star formation (MD99), possibly 
with stars as young a s 100 Myr, and the re is HI gas associ- 
ated with the system (IGallart et al.ll200lT ). This means that, 
unlike the similarly luminous Sculptor, this system was able 
to retain its gas for the whole of its lifetime, be it because 
of a larger potential well and/or a less harsh environment, 
with less stripping and a larger possibility of gas infall with 
respect to a satellite of the MW. 

The most recent study of Phoenix SFH (H09) is based 
on two HST pointings that, even if covering a small fraction 
of the galaxy, reach out to about 4' from the center, giving a 
view of how the SFH changed across the object. The central 
regions contain stars over a range of ages (from 100 Myr 
to the oldest stars), while already at 450pc from the center 
95% of the stars are more than 8 Gyr old. These properties 
appear to be consistent with the star formation region slowly 
shrinking with time (H09). 

In this paper we investigate the global properties of 
Phoenix based on wide-area imaging data collected at the 
ESO VLT with FORS. The observations cover a large re- 
gion, 26'x26', and reach below the horizontal branch of the 
galaxy, containing stars in evolutionary stages representative 
of the whole age range displayed by Phoenix. This allows us 
not only to explore the overall structure of the galaxy but 
also to directly derive the spatial distribution of stars in 
different age ranges. The paper is organised as follows: in 
Sect. [5] we present the observations and the data reduction 
procedure adopted; in Sect. [3] we analyze the overall struc- 
ture of the galaxy; in Sect. [4] we use the tip of the RGB 
to derive an estimate of the distance modulus to Phoenix; 
Sect. [5] deals with the stellar population mix and its spatial 
variations; we conclude with a discussion in Sect. [6] and a 
summary of the results in Sect. [7] 

In a subsequent paper we will complement this pho- 
tometric study with CaT spectroscopy to investigate the 
metallicity and stellar population gradients in this galaxy. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were carried out between June 20 and 
August 20, 2009, in service mode with the ESO VLT in- 
strument FORS2 at UT1 (Antu). FORS2 is a focal re- 
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Parameter 


value 


reference 


notes 


I (y t n n nn . O t o f ^ r ^ f ^ ) 
V*-*-J ZVUU ) u JzUUU/ 


01 h 51 m 06 s , -44°26'42" 


1 




P.A. 


8°± 4° 


2 


beyond 4.5' 




o 3+0 03 


2 


beyond 4.5' 


Rcorc 


1.79' ±0.04' 


2 




Rtidal 


10.56'±0.15' 


2 




Rl/2 


2.30'±0.07' 


2 




(m-M)o 


23.06±0.12 


2 




Distance 


409±23 kpc 


2 




Ly 


0.9 x 10 6 Lq 


1 


referred to a distance of 445 kpc 


V HB 


23.9 


3 




E(B-V) 


0.016 


4 




A v 


0.050 mag 






Aj 


0.024 mag 






1 arcmin/kpc 


0.119 




assuming a distance of 409 kpc 



Table 1. The various rows are from top to bottom: coordinates of the optical centre; position angle P.A., defined as the angle be- 
tween the north direction and the major axis of the galaxy measured counter-clockwise; ellipticity, defined as e = 1 — b/a; King core 
and tidal radius; half-light radius, defined as R1/2 ~ 1-68 X R c , where R e = 1.37±0.01' is the exponential radius; distance mod- 
ulus and heliocentric distance; luminosity in V-band; V magnitude of the horizontal branch; reddening; extinction in Vand /band. 
The last row quotes the conversion factor from arcmin t o kpc at the distance of P hoenix. The P.A . and ellipti city refer to the oute r 
population, aligned north-south. References: 1 = iMateol (l99Sh : 2 = this work; 3 = Holtzman et al ] fcOOCh : 4 = ISchleeel et al ] Gill), 
http:/ /irsa. ipac.caltech.edu/applications/DUS T/ , average value over 2 degrees. The Ay value assumes Ay/E(B-V) =3.1; A7 was derived 
using the ratio A//Ay- = 0.479 from Cardclli et al. ( 1989); note that previous works such as H09 assume larger values of the reddening, 
0.02. 



ducer multi mode instrument that can be used for optical 
imaging, polarimetry, longslit an d multi-object spectroscopy 
|Appenzeller fc Rupprechtlll992t ). It has field of view of 6'8 
x 6.'8 with pixel size (y/25. The detectors are two 2k x 4fc 
MIT CCDs (hereafter referred to as CI and C2). We used 
FORS2 in imaging mode to make a mosaic of 4x4 fields 
covering an area of approximately 26 ' x 26 'centered on the 
coordinates of the optical centre of Phoenix (see Table [1] for 
the galaxy parameters) . One additional field was observed at 
a displaced location from the galaxy as a check on the deter- 
mination of the foreground/background density (see Fig.[T]). 
For each pointing we took 3 dithered exposures of 120s each 
in uhighQ an d 5 dithered exposures of 90s each in Ibcss fil- 
ter (in the subsequent text for brevity we denote filters as 
V instead of vhigh and I instead of Ibcss)- The offset sizes 
between the exposures were such that we covered the gaps 
between the two detectors. Table [2] lists the observation log 
of all scientific fields. 

As a part of the standard calibration plan for FORS2, 
observations of standard star fields in these filters were also 
carried out by ESO service mode observers. The log of all 
standard star observations is in Table [3] 

The basic data reduction steps, consisting of bias sub- 
traction and division by the normalized twilight flat field 
exposure, were performed with the ESO FORS imaging 
pipeline. 

For each pointing such pre-reduced images were first 
aligned to the coordinates of one of the individual expo- 
sures, which was acting as the reference image. The regions 
of the CCD that had no data were masked out - in par- 
ticular the so-called slave CCD (bottom CCD) has a large 



1 This filter has a central wavelength of 555nm and 
a FWHM of 123.2; for the transmission curve see 
http: / / www.cso.org/sci / facilities / par anal / instruments / fors / doc /VLT-MAN-ESO-1 
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Figure 1. Schematic representation of the coverage of our 
VLT/FORS photometric data for the Phoenix dwarf galaxy. The 
black squares indicate the individual pointings, showing the re- 
gions of overlap. The inner dotted ellipse (with parameters e = 
0.4, P.A.= 95°and semi-major axis = 115") denotes the region 
where the elliptical isople ths show a sharp rota t ion in their major 
axes in the analysis of lMarti'nez-Delgado et al.l (l999), and where 
most of the young stars are found. The outer ellipses indicate two 
measurements for the extent of Phoenix: the dotted one (with pa- 
rameters e = 0.3, P.A.= 5°and semi-maj or axis = 15.8') shows the 
nomin al tidal radius from the study of [Martmcz-Dclgad o et al.l 
l|l999l l. and the dashed one (with parameters e = 0.3, P.A.= 
5°and semi-major axis = 8.7') where the ex ponential fit to the 
surfac e brightness profile approaches zero from lvan de Rvdt et al.1 
lll99ll l. 
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Field name 


RA(deg) 


DEC(dcg) 


Date and UT of start of observation 


Filter 


cxptimc [sec] 


airmass 


seeing [arcscc] 


pnoenix-imgOl 


97 ooftO/i 9 
z * ,yyoU4z 


AA 9Qflfi1 


90n0 f!7 97 G-QQ-1/1 f!07 

zuuy-u 1 -z ( y.oo. 14. uy < 


V 


q \/ 1 on 

O X 1ZU 


1 


.073 


n fin 
U.DU 








Ci- A 1 -HQ QQO 

y. 41. Uo.ooz 


1 


k \y nn 
ox yu 


I 


.069 


n fin 
U.DU 


pnocnix-imgOz 


97 Sztfi^JiQ 
A 1 .o4oOoo 


AA 9fififi1 


9000 07 97 in-nn-9Q ^ftQ 
zuuy-u / - z / 1u.uu.zy.000 


v 

V 


Q v 1 90 
O X 1ZU 


1 


.062 


u. Do 








1 O-nfi-9/1 /ICQ 
lU.Uo.Z4.4oo 


1 


^/ on 
o x yu 


I 


.062 


U.oU 


pnoenix-imgUo 


97 701 A^H 


AA 9fifiQQ 


9000 07 98. R^^-^R Q99 

zuuy-u { -zo o.zo.oo.ozz 


v 

V 


o v 1 90 
O X 1ZU 


I 


.138 


U. / U 








o.oo\oy.4oy 


1 

1 


\/ on 
o x yu 




126 


u.oo 


• • 

pnocnix-img04 


97 C^c:^ftft7 
Z ( .00000 ( 


AA 9fi7^ft 
-44. Zo i 00 


9000 fl7 9S S-/1ft-m f^90 

zuuy-Ui-zo O.40.UI.0ZU 


V 


Q v 1 90 
O X 1ZU 


X 


.109 


u.o / 








o.Oo.4y.ooO 


1 

1 


^ ss on 
o x yu 


I 


.100 


U.oo 


pnoenix-img05 


z 1 .yyozo 


AA QGQfiG 

-44.oyooy 


9000 fl7 OQ O-O^-QC; Q^/1 

zuuy-Ui-zo y.uo.oo.o04 


V 


Q v 1 90 

o X IzU 


I 


.091 


U. / U 








Q- 1 Q-9Q 01Q 

y.io.zo.oiy 


1 


v on 
o x yu 


I 


.084 


U.D1 


pnocnix-img06 


97 fi/1 QCQQ 


A A QGQG 

-44.oyoy 


9000 fl7 9/I Q.OQ.Ifl 7QQ 

zUUy-U/-z4 o.Zo.lU./oo 


V 


Q v 1 9n 
O X IZU 


I 


.163 


U. / o 








Q-Qft-nQ ftft7 
O.oO.Uo.OO I 


1 

1 


\/ on 
o x yu 


1 


.149 


U.oo 


• • 

pnoenix-img07 


97 7nneQQ 


A A QQQQc; 

-44.oyoyo 


9000 07 9/i e-nQ-ne nfti 
ZUUy-U(-z4 o.Uo.Uo.UOl 


V 


q ^ i on 
o X IzU 


1 


.216 


U.oo 








fi-1 E^-/10 71 Q 
O.10.4U. ( lo 


1 


^ s/ on 
o x yu 


I 


.187 


U.oo 


pnocnix-imgOo 


97 ^aosq 

Z / .004UOO 


A A QQQQ/1 
-44.03004 


900Q 07 9fi Q-QO-^ft 0.^0 

zuuy-Ui-zo y.ou.oo.oou 


v 

V 


o v 1 90 
O X IZU 


I 


.071 


u.oo 








0-QQ-/1/1 7flt^ 
y.oo. 44. / oO 


1 


^/ on 
o x yu 


I 


.068 


u.oy 


pnoenix-img09 


97 GQft7G9 

z ( .yyo / yz 


-44.0UUOO 


9000 fl7 9fi O-f^l -99 iRQfi 

zuuy-Ui-zo y.oi.zz.ooo 


V 


q ^ i on 
o X IzU 


I 


.065 


u.oo 








y.oy. iz. lor 


I 

1 


^ v 00 

o x yu 


I 


.064 


n fin 

U.OU 


pnocnix-imglO 


97 fi/i 01 9c; 
z / .o4y izo 


-44.0U1 1 1 


90n0 fl7 9/1 fi-/lfi-Q7 f^Qfi 
ZUUy-Ui-Z4 o.4o.Oi.0oo 


V 


q v i on 

O X IZU 


1 


.130 


n fifi 
u.oo 








S-F^ft-/11 Ofl/1 
O.OD.41.yo4 


1 

1 


^ s/ on 
o x yu 


1 


.119 


n fin 
U.OU 


pnoenix-imgl 1 


97 RQfi7^ 

z ( .oyo ( 


a a e;ni 1 ft 

-44.001 IO 


900Q Oft 91 Q-Q^-^/l 111 

zuuy-uo-zi y.oo.04.iii 


V 

V 


Q v 1 90 
O X IZU 


1 


.324 


n fi*3 
u.oo 








y.4o.4o.y4z 


I 

1 


^ v QO 

o x yu 


1 


.298 


n fi^ 
u.oo 


pnocnix-imglz 


97 ^91 ft7 
Z i .OOZ±D / 


a a c;ooqi 

-44.0UU01 


900Q 07 9Q Q-lft-0 1 ^ 79ft 

zuuy-u/-zy y.io.uo./zo 


V 


o v 1 00 
O X IZU 


1 


.079 


n 

U.oo 








9:23:58.611 


1 


5x 90 


1 


.074 


0.84 


phoenix-imgl3 


27.996333 


-44.60359 


2009-07-29 9:35:54.979 


V 


3x120 


1. 


.070 


0.73 








9:43:42.634 


1 


5x 90 


1 


.067 


0.69 


phocnix-imgl4 


27.847792 


-44.6043 


2009-07-29 9:56:11.015 


V 


3x120 


1 


.064 


0.83 








10:04:22.331 


1 


5x 90 


1 


.064 


0.80 


phoenix-imgl5 


27.698542 


-44.60388 


2009-08-20 6:24:45.464 


V 


3x120 


1. 


.199 


0.83 








6:32:43.001 


1 


5x 90 


1 


.182 


0.75 


phocnix-imgl6 


27.551292 


-44.60303 


2009-08-20 6:45:15.295 


V 


3x120 


1 


.157 


0.86 








6:53:06.891 


1 


5x 90 


1 


.143 


0.80 


phoenix-imgl7 


28.246042 


-44.28881 


2009-08-20 7:05:48.935 


V 


3x120 


1. 


.126 


0.84 








7:13:51.822 


1 


5x 90 


1 


.115 


0.64 



Table 2. Table of FORS imaging observation for the Phoenix dT. The seeing is the average stellar FWHM from the final, combined, 
image. 



area that gets no light - and the aligned images were then 
median combined using IRAF task imcombine. The image 
quality was measured on combined images and the FWHM 
values are given in Table [2] We checked that the FWHM 
measured on the combined images was similar to those in 
the individual exposures. 

Since the internal regions of Phoenix are relatively 
crowded, we decided to perform PSF fitting photometry 
on all of our pointings to get uniform photometry over the 
whole area. This was done using the stan dalone version of 
the DAOPHOT II package (jStetsonll 19871 ). 

The routines FIND + PHOTOMETRY + PICK + PSF 
+ ALLSTAR were run separately for each photometric band 
and for the two FORS CCD chips, accounting for the slightly 
different read-out noise of the CCD chips. 

Default parameters were used for the FIND routine to 
deal with the rejection of bad pixels and elongated objects 
along the direction of rows and columns. First, objects were 
located by imposing a 3<r and 4<r threshold above the back- 
ground on the individual images for the V and I band, re- 
spectively. This threshold value was chosen after inspecting 
the number of detected objects as a function of threshold pa- 
rameter, and selecting the value corresponding to the "knee" 



where the number of detected objects is changing rapidly 
with small changes in a for threshold. 

The choice of the PSF stars was done interactively for 
each field, separately for each filter and each CCD. First 
the routine PICK was run to automatically select a sample 
of unsaturated, relatively isolated stars per frame; a PSF 
was then created, and an image containing the residuals us- 
ing ALLSTAR. All the PSF stars were then visually ex- 
amined both on the original image and on the image with 
subtracted stars. Galaxies, blends or stars with close com- 
panions were manually removed from the sample, and the 
PSF re-determined using this subsample of stars. The proce- 
dure was repeated until the only remaining stars to be used 
for the PSF determination were non-saturated and relatively 
isolated stars. Typically, between 15 and 30 such stars per 
frame were used to determine the PSF. We used the op- 
tion for a spatially non variable PSF, after testing that no 
significant differences were introduced by allowing for a spa- 
tially variable PSF. The analytical first approximation used 
as model for the PSF was a Moffat function. 

The above procedure, except the creation of the PSFs, 
was then re-run on the residual images, in order to extract 
stars that were not previously found due to the wings of 
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Field name Date and UT of observation Filter airmass exptime [sec] seeing [arcsec] 



E5 


2009- 


■06- 


20 


22: 


42: 


09. 


728 


V 


1 


.071 


3 





.53 










22: 


43: 


39 


.426 


I 


1 


.071 


1 





.45 


E7 


2009- 


•07- 


24 


02: 


02: 


19. 


.507 


V 


1 


.067 


3 





.80 










02: 


03: 


50. 


.096 


1 


1 


.067 


1 


1 


.00 


L92 


2009- 


•07- 


27 


10: 


27: 


10. 


238 


V 


1 


.163 


3 





.50 










10: 


29: 


05 


.406 


1 


1 


.165 


1 





.53 


MarkA 


2009- 


■07- 


■27 


05: 


:58; 


03. 


331 


V 


1 


.055 


3 





,15 










05: 


59: 


.28. 


.549 


1 


1 


.056 


1 





,13 


PG0231.1 


2009- 


•07- 


■24 


10: 


23: 


00. 


265 


V 


1 


.176 


3 





.65 










10: 


24: 


:28. 


.883 


1 


1 


.175 


1 





.75 


PG0231.2 


2009- 


■07- 


28 


10: 


13: 


18. 


350 


V 


1 


.171 


3 





.60 










10: 


14: 


13 


.538 


1 


1 


.169 


1 





.65 


PG1525 


2009- 


•06- 


■20 


22: 


50: 


50. 


.644 


V 


1 


.597 


3 





.55 










22: 


52: 


19. 


.402 


1 


1 


.585 


1 





.50 


PG1657 


2009- 


■07- 


■27 


00: 


32: 


21. 


166 


V 


1 


.213 


3 





,53 










00: 


33: 


46. 


.404 


I 


1 


.212 


1 





,13 


PG2213.1 


2009- 


•07- 


■27 


07: 


10: 


.52. 


208 


V 


1 


.108 


3 





.50 










07: 


12: 


.17. 


.476 


1 


1 


.109 


1 





.18 


PG2213.2 


2009- 


■07- 


■28 


10: 


22: 


26. 


352 


V 


2 


.015 


3 





.73 










10: 


23: 


51 


.520 


1 


2 


.035 


1 





.65 


PG2213.3 


2009- 


•07- 


■29 


06: 


17: 


34. 


.819 


V 


1 


.100 


3 





.50 










06: 


49: 


00. 


.078 


1 


1 


.101 


1 





.15 


PG2213.4 


2009- 


■07- 


29 


09: 


11: 


13. 


858 


V 


1 


.437 


3 





.78 










09: 


12: 


39 


.046 


1 


1 


.445 


1 


1 


.00 


PG2213.5 


2009- 


•07- 


29 


10: 


29: 


37. 


.794 


V 


2 


.181 


3 


1 


.50 










10: 


31: 


03. 


.072 


1 


2 


.204 


1 





.88 


PG2213.6 


2009- 


■08- 


20 


07: 


28: 


30. 


579 


V 


1 


.359 


3 





.78 










07: 


29: 


09 


.972 


I 


1 


.362 


1 





.70 


T.Phe-1 


2009- 


•07- 


28 


09: 


27: 


24. 


.530 


V 


1 


.088 


3 





.60 










09: 


:28: 


49. 


.538 


1 


1 


.089 


1 





.13 


T.Phe.2 


2009- 


•07- 


29 


10: 


23: 


12. 


828 


V 


1 


.149 


3 





.75 










10: 


24: 


38. 


.096 


I 


1 


.151 


1 


1 


.38 


T.Phe.3 


2009- 


•08- 


20 


07: 


33: 


51. 


860 


V 


1 


.079 


3 





.68 










07: 


34: 


30 


.794 


I 


1 


.079 


1 





.60 



Table 3. Table of FORS imaging observation for standard fiels. The seeing is the average stellar FWHM from the final, combined, image. 



brighter stars. We checked that this procedure was not fit- 
ting noise by overplotting the position of these fainter stars 
on top of the image of residuals. 

Finally, DAOMATCH and DAOMASTER were used to 
combine V and I ALLSTAR PSF fitting photometry of each 
pointing - separately for the 2 CCD chips. Only objects de- 
tected in both bands and with magnitude errors < 0.5 mag 
were retained. 

Aperture corrections were computed for each field, chip 
and filter separately using aperture photometry and a curve- 
of-growth method on a set of well exposed and isolated PSF 
stars. 



2.1 Astrometry and Photometric calibration 

The J2000 celestial coordinates of the detected objects 
were derived with CataXcori0. The astrometric solution was 



2 This code was developed by Montegriffo at INAF- Osservatorio 
Astronomico di Bologna, with the aim of cross-correlating cata- 
logues and finding astrometric solution. The code is routinely used 
in publications f rom staff of the INAF - Osservatorio Astronomico 
di Bologna (e.g. iBellazzini et alj|201ll) . 



found by fitting typically a polynomial of third degree be- 
tween at least 46 stars per ponting and the reference cata- 
logue (USNOA2 and GSC2). The average r.m.s. of the solu- 
tion was 0.2"in both RA and DEC. 

Given that the standard fields used are not crowded, we 
performed aperture photometry rather than PSF photome- 
try. We located only luminous standard stars, by adopting a 
threshold of 15a above the background on the individual im- 
ages 0. The aperture at which to calculate the instrumental 
magnitude was derived by examining the run of the magni- 
tude as a function of aperture size for those bright stars that 
are neither saturated, nor close to the border of the image, 
and for which the curve of growth reaches a constant value. 

In general an aperture of 25 pixels radius was found to 
include all the stars' flux for most of the fields of standard 
stars. However, this aperture radius was too large for a rela- 
tively crowded field such as E7; furthermore only a handful 
of stars passed the above criteria for fields E5 and PG0231: 
in these three cases we therefore used a smaller aperture 



3 We checked that reducing this threshold would not significantly 
increase the number of detected stars with catalogued magnitude 
per frame nor the quality of our calibration. 
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Night & CCD chip ZP V-band ZP I-band 



21 


Jun CI 


-27 


.934+ 





.025 


-27.422+ 


0. 


.043 


21 


Jun C2 


-27. 


.997+ 


0. 


009 


-27.400± 


0. 


044 


21 


Jul CI 


-28. 


.035+ 


0. 


.007 


-27.418± 


0. 


007 


21 


Jul C2 


-28. 


.022+ 


0. 


.007 


-27.410± 


0. 


Oil 


27 


Jul CI 


-27. 


.956+ 


0. 


038 


-27.485± 


0. 


042 


27 


Jul C2 


-28. 


.027+ 


0. 


.010 


-27.394± 


0. 


008 


28 


Jul CI 


-28. 


.012± 


0. 


030 


-27.415± 


0. 


.027 


28 


Jul C2 


-28. 


.049+ 


0. 


017 


-27.373± 


0. 


023 


29 


Jul CI 


-28. 


.080+ 


0. 


.025 


-27.407+ 


0. 


.019 


29 


Jul C2 


-28. 


.032+ 


0. 


026 


-27.393+ 


0. 


033 


20 


Aug CI 


-28. 


.005+ 


0. 


034 


-27.378± 


0. 


025 


20 


Aug C2 


-28. 


.049+ 


0. 


049 


-27.439± 


0. 


051 



Table 4. Zero-points (and their errors) of the photometric cali- 
bration transformations (Eq. [1} derived for each night of obser- 
vations, separated per band and CCD chip. 



radius and derived the aperture correction choosing appro- 
priate isolated stars manually. 

Applying the above described procedure, we extracted 
13 standard stars for the night of June 21 (9 chipl + 4 
chip2); 91 (61 + 30), 58 (36+22), 19 (11+8), 28 (19+9) for 
July 24,27,28,29 respectively; 22 (12+10) for August 20. 

We note that the data available for the standard stars 
are not sufficient to fit simultaneously the zero point, the ex- 
tinction coefficient and the color term for each band, night 
and chi p separately. We there fore follow an iterative ap- 
proach l|jerien fc Re ikuba 2001), which we apply to the two 
CCD chips separately. 

For each band, we use the equation 

m inst r = M + Z m + k m x X m + c m x (V - I) (1) 

where mi nstr and M are the instrumental and calibrated 
magnitude, X m the airmass of the pointing, Z m the zero 
point, k m the extinction coefficient and c m the color term. 
The calibrated photometry for the s tandard stars was ex- 
tracted from the online catalogues of IStetsonl ([2000 

On a first pass, we derive the zero point and extinc- 
tion coefficient for each night separately by fitting the in- 
strumental magnitudes and keeping the colour term fixed 
to zero. The derived coefficients are then applied to the 
calibrated magnitudes of the respective nights, obtaining 
m stcp i = M + Z m ,fl t + fcm.fit x X m ; the data of all the nights 
are then joined together to increase the statistics, and a color 
term is fitted to mi nstr — m stcp i = c m ,fit x (V — I) + const. 
The fitting procedure of the zero point and extinction term 
night by night is then repeated, this time by keeping the 
color term fixed to Cm, fit- Any additional dependency on 
the color has been added to the color term. 

The average color term in V-band is -0.063+0.012 for 
CI and -0.053+0.006 for C2, and is consistent with zero for 
the I-band, 0.010+0.006 for CI and 0.001+0.01 for C2. This 
compares well with the color terms available on the ESO 
quality control webpages for P8J0, when available for the 
color used here. The average extinction term in V-band is 

4 http:/ /www3.cadc-ccda.hia-iha.nrc- 
cnrc.gc.ca/community/STETSON/standards/ 

5 http:/ /www. eso.org/observing/dfo/quality/FORS2/qc/photcocff/ 



0.141+0.015 for CI and 0.134+0.008 for C2, while in I-band 
it is 0.055+0.011 for CI and 0.039+0.010 for C2; the deter- 
minations for the two chips are consistent with each others 
within la. The nightly zero points are listed in Table(4] The 
calibrated magnitudes for the science targets were finally de- 
rived by iteratively applying the calibration transformations 
to the instrumental magnitudes, to which aperture correc- 
tions had been previously applied. The errors on the de- 
rived aperture corrections and coefficients of the calibration 
transformation are included in the error on the calibrated 
magnitudes. 

We used the region of overlap between adjacent fields 
to place our photometry on a common photometric scale: 
for this we applied the magnitude shifts derived by com- 
paring the magnitudes of objects in the overlap regions to 
the various pointings and tied all of them to the system of 
pointing 06. The applied shifts were typically < ±0.1mag; 
in part the magnitude shifts applied are the consequence 
of the errors in the calibration of the standard fields, but 
mostly they arise from t he imperfect fiat-fieldi ng due to sky 
concentration in FORS (|Freudling et al.ll2007l n The errors 
in the determinations of the magnitude shifts were propa- 
gated onto the magnitude errors of the individual objects. 
Obviously pointing 17, which has no overlap with the main 
mosaic region, could not be placed on the same photometric 
system. Since however we use this pointing only for checks 
on the determination of the foreground/background density, 
this will not affect our analysis. 

The final step is to merge the catalogues for the various 
pointings in one catalogue containing a unique set of ob- 
jects. This is done by performing weighted averages on the 
observed quantities for stars with multiple measurements. 

At this stage we also evaluate which objects we consider 
as "stars" on the basis of their PSF fitting shape parame- 
ters and errors in magnitudes. Figure [2] shows the run of 
the DAOPHOT photometric parameters "sharpness" and 
"chi", as well as the overall magnitude error as a function 
of the calibrated magnitudes for the sample of unique (com- 
bined) measurements. We exclude a priori those objects with 
relatively large magnitude errors, i.e. av > 0.3 mag and 
oi > 0.4 mag. We then examine the distribution of sharp- 
ness versus magnitude for this sub-sample of stars, deriving 
its scatter per bin of magnitudes. Given the asymmetric dis- 
tribution around null values of the sharpness parameter, we 
retain as bona-fide stars those objects that fall within the 
region defined by the exponential envelopes fitted to +3x 
and — 2x the scatteiQ. The catalogue of bona-fide stars con- 
sists of 19286 stars including the control field and 18919 in 
the main mosaic area (see Table [5]). 

Treating photon counts of each star as a Poisson vari- 
able, a S/N=10 would correspond to an error of magnitude 
of ~0.1mag, and a S/N =5 to ~0.2mag. This would imply 

photcocffs_fors2.html; note that, as stated on this webpage, the 
errors on the tabulated quantities are underestimated. 

6 Although the report from Freudling et al. described in detail 
the photometric accuracy of FORS1, it is expected that FORS2 
flat-fields present similar effects as the two cameras have the same 
optics and differ only in the detector type 

7 The scatter is derived in an iterative manner and it is the scaled 
median absolute deviation (m.a.d.) of the sharpness per magni- 
tude bin of 0.5mag. 



Survey of the Phoenix transition type 7 



V) 
V) 

<v 
c 

Q. 
i_ 
D 
_C 



1 - 



r 



-1 



6 - 



x 



2 - 







0.5 r 



0.4 t 



sr 0.3 



0.2 r 



0.1 r 



0.0 



- 1 — i — i — i — i — i — i — 



i — 1 — 1 — 1 — i — 1 — • — 1 — r 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 




H 1 I— 1 1 1 H 



H 1 1 1 1 1 1 1 1 1 1 1 1— +- 




-H-<- 





H 1 1 1 H 



H 1 1 1 1 1 1 1 1 1 1 H 





16 18 20 22 24 26 16 18 20 22 24 26 

V I 

Figure 2. Photometric parameters for the detected sources (cyan points) as a function of the calibrated V- (left) and I-(right) magnitudes. 
The top and central panels show the sharpness and \ parameters, respectively, from DAOPHOT. The bottom panels show the photometric 
errors. Beside the magnitude errors output from DAOPHOT, these errors factor in the uncertainties from PSF correction, photometric 
calibration and the shifts applied to the individual magnitues so as to bring all the photometry to the system of pointing 6. The black 
points show only those sources which we consider as "stars" on the basis of their sharpness parameter, and that have ay < 0.3 mag and 
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Name 


RA [deg] 


DEC [dcg] 


V 


cry 


I 


07 


X 


sharpness 


phxOOOOl 


27.9264897 


-44.2962170 


25.0495 


0.1396 


23.8995 


0.1396 


1.2195 


0.6250 


phx00002 


28.0484214 


-44.2960092 


25.4608 


0.1679 


23.8604 


0.1679 


1.1445 


0.6710 


phx00003 


28.0380659 


-44.2959729 


25.3526 


0.0983 


22.8875 


0.0983 


0.7605 


0.2435 


phx00004 


28.0376092 


-44.2955426 


24.8144 


0.0926 


24.5143 


0.0926 


0.6230 


0.2570 


phx00005 


28.0383987 


-44.2954079 


25.3673 


0.1125 


23.0559 


0.1125 


0.7450 


0.1105 


phx00006 


27.9966041 


-44.2952897 


24.2082 


0.0921 


23.2352 


0.0921 


1.4590 


0.5680 


phx00007 


27.9471308 


-44.2953097 


24.5857 


0.0768 


24.0239 


0.0768 


0.7310 


0.1025 


phx00008 


28.0220882 


-44.2952300 


25.6432 


0.1224 


23.3985 


0.1224 


0.9160 


0.4125 


phx00009 


27.9638667 


-44.2951593 


24.3899 


0.0847 


23.6765 


0.0847 


0.9075 


0.3655 


phxOOOlO 


28.0697991 


-44.2950072 


24.5369 


0.0996 


23.6116 


0.0996 


1.2930 


0.5940 



Table 5. Catalogue of bona-fide stars from FORS photometry of the Phoenix dT, including pointing 17. The columns list the name of 
the star (1), RA (2) and DEC (3) in degrees, V magnitude (4), error in V magnitude (5), I magnitude (6), error in I magnitude (7), x 
(8) and sharpness (9) parameters (8). This is a sample of the table; the complete version will be published online. 
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Figure 3. Magnitude difference for stars in the overlapping regions between different pointings and with double measurements as 
a function of calibrated magnitudes (left: V band; right: / band). Here we plot only the stars that passed the selection criteria in 
sharpness and photometric errors (see Fig.[2j- The scaled median absolute deviation at the bright end side of the distribution is 0.054 
mag for the V-band (19 < V < 21) and 0.065 mag for the I-band (18 < I < 20), corresponding to an error in the individual measurement 
of 0.038 mag and 0.046 mag. At the faint end side these numbers are 0.14 and 0.13 (24 < V < 25 and 23 < I < 24. respectively), 
corresponding to error in the individual measurement of 0.1 and 0.09 mag. At the bright end the photometric errors are 0.033 and 0.035 
for Vand /band, respectively, and at the faint end they are 0.076 and 0.069 at the faint end. Therefore, the magnitude differences from 
stars with multiple measurements are inflated by an extra 0.01 and 0.03 mag, and at the faint end by 0.06 and 0.06 mag. 



for this data-set a S/N=10 at the magnitude Vio ~ 24.8, 
Jio ~ 23.6, correspondig to 11737 stars in the main mo- 
saic area. The magnitude limits at S/N =5 are Vg ~ 25.4, 
Is ~ 24.0. Here we have quoted the shallowest of the val- 
ues among all fields and chips. The number of stars in the 
main mosaic area above the S/N=5 and S/N=10 magnitude 
limits is 15238 and 11737, respectively. 

Figure shows the magnitude differences of the bona- 
fide stars with double measurements. These differences will 
be used to estimate the typical magnitude and color error 
in various magnitude bins (see Sect. 15.1]) : in these values are 
therefore included internal photometric errors, errors from 
aperture corrections, errors due to the photometric calibra- 
tion applied to different pointings, and effects of crowding 
(including possible mismatches of stars in the most crowded, 
inner regions). 

In order to understand in more detail the role of crowd- 
ing in the error budget and on the completeness limit from 



the inner to the outer parts, we also perform artificial star 
tests on pointing 06. This is representative of the four cen- 
tral, most crowded pointings (see Fig. [T]). Since the FORS 
f.o.v. is 6.8' x 6.8', one single pointing already extents from 
the center of Phoenix to about half nominal tidal radius 
(VDK91, MD99). 

We add artificial stars to the combined V and I images 
for CI and C2 separately, using the PSF-model derived for 
pointing 06. The coordinates and magnitudes are randomly 
chosen from a uniform distribution, and cover the whole ex- 
tent of the pointing and instrumental magnitudes of stars 
in Phoenix (corresponding to calibrated magnitude in the 
ranges F[20,27] and /[18.5,26.5]). We inject 200 artificial 
stars per band per chip (about 5-10% of the observed num- 
ber of stars in these chips), and repeat the experiment 5 
times. Finally we apply the same reduction and selection 
procedure that we used to produce the final catalogue of 
stars. We find that crowding does not significantly affect 
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Figure 4. Hess diagram in bins of 0.3'of the spatial distribution 
of the stars down to V = 24.8, I = 23.6 (S/N> 10). The solid 
ellipse shows the nominal tidal radius derived in this work, rt = 
10.56'±0.15' (e = 0.3, P.A.= 5°); for the inner dash-dotted ellipse 
see caption of Fig. [l] 

the 50% completeness limit (i.e. the magnitude at which the 
fraction of recovered versus injected number of stars is 0.5): 
in /band this is found to be constant around J~23.5-23.8 
over the explored range of R, while in Vis V~24.2 in the 
inner 1.5', and at larger R it remains constant at V~~25. 
These values are very similar to the S/N= 10 limits quoted 
above. 

The difference of recovered and input magnitudes of the 
artificial stars as a function of recovered magnitude is fairly 
symmetric around zero, showing that - except for a handful 
of stars - blending is not affecting the recovered magnitudes. 
The scaled m.a.d. between the difference of input and recov- 
ered calibrated magnitude is smaller than those determined 
using the stars from double measurements, e.g. in the most 
crowded region is 0.07 mag for ^between 24,25 and 0.09 mag 
for ^between 23,24. The errors from Fig. [3] which we use 
throughtout the analysis, can then be considered as conser- 
vative. 

3 THE OVERALL STRUCTURE OF PHOENIX 
3.1 2D distribution 

When analyzing the 2D distribution of stars in Phoenix, 
MD99 found a sharp variation in the trend of the ellipticity 
(e = 1 — b/a, where b and a are the minor and major axes of 
the ellipse) and position angle (P. A.) of the isophleths with 
radius, with e changing from 0.4 to 0.3 and the PA. from 
95° to 5° around 115". The authors noted that the majority 
of young stars were contained in the flattened central "com- 
ponent" tilted with respect to the main body of the galaxy 
and the authors proposed this may be a disk-halo structure, 
as for large spiral galaxies. 



Figure 5. The variation with projected radius of the central po- 
sition (x and y coordinate, a) and b), respectively), c) ellipticity 
and d) position angle for the entire stellar population of Phoenix, 
for the stars brighter than V = 24.8, I = 23.6 (S/N> 10) from our 
FORS imaging data. The variation of the central position is with 
respect to the value listed in Mateo (1998). The horizontal lines 
in panel c) and d) indicate the values of the ellipticity and P. A. 
found for the outer Phoenix population by Martinez-Delgado et 
al. 1999. 



We construct Hess diagrams of the spatial distribu- 
tion of Phoenix stars, using the sample at S/N > 5 and 
10 (brighter than V 5 ~ 25.4, 7 5 ~ 24.0, and than Vio ~ 24.8, 
iio ~ 23.6, respectively), in bins of 0.3'. Since in some parts 
there appear to be voids or enhancements in the stellar num- 
ber counts with respect to the surroundings, e.g. in the over- 
lap region between pointings, we prefer not to smooth the 
spatial distribution because this may enhance such features. 
Figure [4] shows the resulting distribution for the S/N>10 
sample: it is clearly visible that the inner parts display a 
different distribution than the outer parts, being almost per- 
pendicular to each other. 

We run the IRAF task ELLIPSE on the above Hess di- 
agrams, letting the center, e and PA. of the isophotes at 
different radii as free parameters. The results are shown in 
Fig.[S]for the S/N>10 sample: within a radius of 1 and 2 ar- 
cmin, the center appears to be located about 0.6 -0.2 arcmin 
east a nd 0.1 arcmin south from the value listed in iMateol 
|P998j); the ellipticity decreases from a value of 0.44±0.05 at 
l'to 0.1 at about 3', to increase again, with an average value 
of 0.28±0.03 at R ^ 4.5'; the PA. varies from about 80° 
within 2'and decreases reaching an average value of 8° ± 4° 
at R ^ 4.5'. These results are very similar to those resulting 
from the S/N> 5 sample. They are also in good agreement 
with those of MD99, even though their results on the 2D 
distribution may have been more affected by crowding be- 
cause of the larger seeing of the observations. We discuss the 
possibility of this being a disk-halo structure in Sect. [5] 

The intermediate values of ellipticity and PA. that can 
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Table 6. Values of the density of contaminants [number 
arcmin -2 ] as derived from the weighted mean of the outer points 
of the surface number profile of the various populations (column 
2), from the objects at elliptical radii >12. 6' (column 3) and from 
pointing 17 (column 4). Note that since the MS display an asym- 
metric distribution we do not derive surface number profiles for 
them. For those populations that have stars as faint as the V5 , 1$ 
magnitude, we perform two determinations: one considering the 
stars brighter than the Vfo, /10 limit (labelled as 10) and those 
brighter than the V5, I5 one (labelled as 5). Note that the major- 
ity of the determinations are compatible within lcr between the 
3 methods. 





weighted mean 


R >12.6' 


field 17 


ms 


X 


0.012 ± 0.006 


0.022 ±0.022 


ms 10 


X 


0.0 


0.022±0.022 


bl 


O.llzt 0.02 


0.092 ± 0.017 


0.087±0.04 


rc 


0.543± 0.055 


0.474 ± 0.038 


0.281±0.078 


agbbump 


0.056±0.018 


0.037 ± 0.011 


0.065±0.037 


rhb 


0.338± 0.046 


0.361 ± 0.033 


0.238±0.072 


bhb 5 


0.07± 0.023 


0.092 ± 0.017 


0.043±0.031 


bhb 10 


0.053± 0.020 


0.049± 0.012 


0.022±0.022 


rgb 


0.055± 0.020 


0.046 ± 0.012 


0.022±0.022 


all 5 


8.44 ± 0.24 


8.15 ± 0.16 


7.40±0.40 


all 10 


5.52 ± 0.19 


5.35±0.13 


5.58±0.35 



be seen around 2' are probably due to the super-position of 
the inner and outer component. 




Elliptical radius [arcmin] 

Figure 6. The surface number density profile for the Phoenix 
dwarf galaxy with overlaid best-fitting King, Sersic, exponen- 
tial and Plummer models (solid, dashed, dotted, dash-dot lines, 
respectively). The density of contaminants, 5.52 ± 0.19 stars 
arcmin" 2 calculated from a weighted average of the outer 3', has 
been subtracted from each point. The error bars are obtained by 
summing in quadrature the errors from Poisson statistics and the 
error in contamination. The central point was excluded from the 
fit. Our best fit to the data is a Sersic profile with Sersic radius, 
Rs = 1.82' ±0.06' and m= 0.83±0.03 with a reduced \ 2 = 1-8; a 
King profile with core radius r c = 1.79' ± 0.04' and tidal radius 
rt = 10.56' ± 0.15' gives a reduced \ 2 = 3.5;. The best-fitting 
parameters for the different models are summarized in Table 171 . 



3.2 Surface number density profile 

There have been previous determinations of the surface 
number density profile of Phoenix, with results rather differ- 
ent from each others: VDK91 used a scanned ESO/SRC II- 
IaL plate and found that the surface brightness profile is not 
well fit by a King profile but decreses exponentially reach- 
ing the zero level at 8.7' from the center; on the other hand, 
MD99 find that the profile can be well fit by a King pro- 
file, and determine a nominal tidal radius of rt = 15.8^2"! 
arcmin. 

We derive the surface number density profile of the over- 
all stellar population of Phoenix using both the data-sets 
with S/N> 5 and 10. We sample the number counts in ellip- 
tical binijf]. For the values of ellipticity of the elliptical bins 
we adopt the values for the outer components of MD99, i.e. 
e = 0.3 and P.A.= 5°, given the good agreement between our 
results and those of MD99, and the fact that MD99 derived 
those parameters from a single plate rather than a mosaic 
of images. 

Given the discrepant measurements of the extent of 
Phoenix in the literature, we decided to estimate empir- 
ically the density of contaminant objects (MW foreground 
stars and unresolved galaxies), rather than relying on uncer- 
tain selection of a region supposedly free of Phoenix stars. 
The following two approaches were adopted: a) we compute 



8 Throughout the manuscript the projected elliptical radius of a 
point (x, y) is defined as R = y/ x 2 + j/ 2 /(l — e) 2 , where e is the 
considered ellipticity, and the galaxy is assumed to be centred on 
the origin, with its major axis aligned with the x-axis. 



the number density profile in elliptical annuli whose semi- 
major axis is <13', i.e. the whole ellipse is contained in our 
main mosaic area; we analyze the outer parts of the profile, 
and adopt as density of contaminants a weighted average of 
the outer points, where the profile becomes approximately 
flat; in this case this is the outermost 3'; b) we analyze the 
number density profile along the projected minor axis of the 
galaxy; this is because in that direction our mosaic covers 
a region corresponding to rather large elliptical radii, 18.6', 
i.e. a region beyond the largest literature measurement of the 
nominal tidal radius of the galaxy (15.8', MD99); along this 
axis, we find that the profile flattens at ~8.8', correspond- 
ing to an elliptical radius of 12.6'; we therefore expect the 
great majority of the objects found at elliptical radii larger 
than 12.6' to be composed by contaminants; we use then the 
objects located at elliptical radius >12.6'to determine the 
density of contaminants. The values we obtain with the two 
methods are consistent with each other (see Table [BJ. As a 
consistency check we compare these determinations to the 
ones derived from pointing 17 (see Table [BJ, i.e. the point- 
ing displaced from the rest of the mosaic area, and also find 
good agreement; however, given the small area covered by 
pointing 17, we will be using the other two methods in the 
following analysis. 

We subtract the density of contaminants derived 
from methods a) and b) from the number density pro- 
file of Phoenix sampled in bins of 0.5' and we fit this 
contamination-subtracted profile to several surface bright- 
ness models using a le ast-square fi t to the data. We used an 
empirical Kin g profile l|King[l96 2l , an exponentia l profile, a 
Sersi c profile ([Sersic 1968) and a Plummer model (|Plummerl 

EaI3). 
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Table 7. Parameters of best-fitting King model (core radius, r c , tidal radius, rt), Sersic model (Sersic radius, ijg, shape parameter, m), 
exponential model (scale radius, rg x h) and Plummer model (scale radius, b) for various stellar populations in Phoenix dSph. For each 
model we list the corresponding reduced x 2 , X 2 e< j = X 2 /(N — v), where N is the number of points in the surface number count profile 
and v the number of degrees of freedom. The errors in the parameters are from Montecarlo simulations but agree well with the formal 
errors from the fitting procedure. These values of the parameters were derived by using the density of contaminants from the weighted 
average of the last points of the surface number profile. The best-fitting profiles have good values of the reduced \ 2 . 
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lar when using the contaminant density derived from either 
method a) or b). In the following we will adopt method 
a) because the fits to the surface brightness profile pro- 
duce slightly better x 2 values; method a) will be replaced 
by method b) only when analyzing the young main se- 
quence stars because of their asymmetric distribution (see 
Sect. [52). 

The results of the fit for the S/N>10 sample are shown 
in Fig. [5] and summarized in Table \7\ In our analysis, the 
Sersic profile is the one that fits the data best, yielding a re- 
duced x 2 = 1-8 for a scale radius R s = 1.82' ± 0.06' and 
m — 0.83 ± 0.03. An exponential profile with exponen- 
tial radius 1.37'±0.0l'(reduced x 2 = 3.2) and a King pro- 
file with core radius r c — 1.79' ± 0.04', tidal radius r t = 
10.56' ± 0.15' (reduced \ 2 ~ 3.5) also reproduce the data 
rather well; the Plummer profile instead clearly overpredicts 
the surface number density at R>7'. The results from the 
S/N>5 sample are very similar, both in terms of shape pa- 
rameters and best-fitting profiles. Because of the intrinsic 
uncertainty in determination of crowding corrections and 
because we are mostly interested in the large scale prop- 
erties of Phoenix, we are not applying a crowding correc- 
tion factor. Instead, we have performed the fit excluding 
the central point of the surface number count profile. It was 
anyway shown in Sect. 12.11 that the effects of crowding are 
not significant. This is confirmed by the results of the fit to 
the surface number count profile derived from the sample of 
stars brighter than the 50% completeness limit at R <1.5': 
the best-fitting parameters are in very good agreement with 
those derived from the sample S/N>10, and they yield very 
similar x 2 values. 

In order to assess whether Sersic and King profiles per- 
form better because of one more free parameter when com- 
pared to exponential an Plummer profiles, we check how the 
performance of the best-fitting profiles would be ran ked us- 
ing the Akaike information criterion (jAkaikeill973l ) in the 
form x 2 + 2 x k, where k is the number of free parameters 
in the fit: the ranking of the best-fitting profiles remains 




The King model has been extensively used to describe 
the surface density profile of dSphs. 

/Kl//i- /„K j 1 - , 1 1 (2) 

It is defined by 3 parameters: a characteristic surface density, 
/o,k, core radius, r c and tidal truncation radius, rt. 

We note that excesses of stars ha ve been found beyond 
the King tidal ra dius in Fornax (e.g. I Coleman et al.l 2005; 
Battaglia et alj I2006I) . Scu lptor (e.g. IColeman et al.l 
20051 : iBattaglia et ail 120081 ) an d in a number o f 



other LG d Sphs, e.g. Cari n a dMaiewski et all l2005h . 
Draco (e.g. IWilkinson et all I2004I). Ursa Minor (e.g. 
Irwin fc Hatzidimitrioul 1 19951 : [Martmez-Delgado et alj 



2001] ). Since the tidal radius is set by the tidal field of 



the host galaxy, the observed excesses of stars have been 
interpreted as tidally stripped stars. However, if dwarf 
galaxies are embedded in massive dark matter haloes, this 
tidal radius loses its meaning of tidal truncation radius 
(we will therefore refer to it as "nominal tidal radius") and 
such excesses of stars need not be tidally stripped stars 
but can simply due to the King model not providing the 
best representation of the surface number count profiles of 
dSphs at large radii. 

The Sersic profile is known to provide a good em- 
pirical formula to fit the projected light distribution of 
elliptical galaxie s and the bulges of spiral galaxies (e.; 



elliptical galaxies and the bulges ol spiral galaxies (e.g . 
Caon et al J 1 19931: ICaldwelllll99a [Graham fc Guzman! 12003! : 



Truiillo et aiT l2004) and also provides a good representa- 
tion of the number surfa ce density of some dSphs (e.g. 
IBattaglia et~ai1l2006l , [2008T ): 



7s (R) = Io.s exp 



R \ 1/m 



(I) 



(3) 



where Jo,s is a scale surface density, Rs is a scale radius and 
m is the surface density profile shape parameter. 

We find that the best-fitting parameters are very simi- 
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the same as according to the reduced x 2 values, therefore 
hereafter we consider only the reduced \ 2 values. 

Our determination of the tidal radius r< = 10.56' ± 
0.15' is smaller than the one from MD99, r t = 15.8^:8 
arcmin; note that the value from MD99 was derived us- 
ing the surface number density between projected radii of 
approximately 2' and 5', and extrapolating the behaviour 
at larger distances, while our determinations comes from a 
much wider area. 

Since the outer parts of the object are almost equally 
well fit by a range of profiles, continuous and truncated, it is 
therefore difficult to establish whether Phoenix was tidally 
truncated by the Milky Way or not. 

In the following we will use the King tidal radius we 
obtained as an indication of the extent of Phoenix, and we 
will refer to it as the "nominal" tidal radius as it is unclear 
whether this is the radius of tidal truncation of the galaxy. 



4 TIP OF THE RGB 

We derive the tip of the RGB of Phoenix in /, Jrgbt, and 
compare it to the determinations from previous works. In 
addition to providing an independent distance determina- 
tion, this also offers a possibility to check our photometric 
calibration. 

For this, we convolve the luminosity function (LF) in 
the I band for stars w ith 18 < I< 21 with a Sobel edge- 
detection kernel [-2,0,2] (|Lee et al.ll 19931 ). which gives a max- 
imum where the discontinuity in the luminosity function is 
greatest (see Fig. [7] left). 

In order to factor in the uncertainties due to the binning 
of the LF, we perform the determination over 100 random 
realizations of the luminosity range used for the calculation 
- that is, we allow the brightest end of the range to vary be- 
tween / 18 and 19, and the faintest end between 21 and 22. 
For each of these random realizations, we perform the convo- 
lution using bin sizes decreasing from 0.15 mag to 0.05 mag 
with step of 0.01 mag. Since there are several discontinuities 
in the LF due to various features in the colour-magnitude 
diagram (see below), we restrict our search for the maxi- 
mum to the region 18 < I< 20; this choice comes from an 
educated guess from a visual inspection of the CMD, and 
from literature values of Jrgbt- 

We find /rgbt = 19.08 ± 0.06, where the error is the 
scaled m.a.d. of the distribution of values. This agrees with 
previous measurements such as those by H09, who find 
/rgbt = 19.14 ± 0.02, and MD99 Irgbt = 19.00 ± 0.07 
mag. These authors use a slightly different version of the 
Sobel filter, i.e. with a kernel [1,2,0,-2,-1], but we checked 
that the use of this kernel makes no difference to our deter- 
mination. 

We adopt a reddening value E(B-V)= 0.016 (see Ta- 
ble [T]), derived from the reddening maps of our Galaxy avail- 
able at http://irsa.ipac.caltech.edu/applications/DUST/ 
centered on the coordinates of the Phoenix dwarf galaxy. 
This corresponds to an ex tinction in V Ay = 0.05 mag 
and in /band Ai = 0.024 jCardelli et al.lll989l ). yielding a 
dereddened magnitude for the RGB Tip /rgbto ~ 19.06. 

As shown in Lee, Freedman & Madore (1993) the abso- 
lute magnitude of the tip of the RGB in I-band changes by 
less than 0.1 mag around the value Mj = —4.0 in the metal- 



licity range -2.2<[Fe/H]< —0.7, which includes the metallic- 
ity range expected for stars in Phoenix (see H09) . Therefore 
by adopting an absolute magnitude Mi^rgbt = —4.0 and 
an extinction Aj = 0.024, a Irgbt = 19.08 would give a dis- 
tance modulus (m-M)o = 23.06±0.12; the error is due to the 
combination of the estimated error on the determination of 
/rgbt and the fact that we are neglecting the (weak) depen- 
dence of the absolute magnitude of the RGB tip in I-band 
with [Fe/H] (i.e. the O.lmag variation of the magnitude of 
the tip around the value Mi = —4.0 mentioned above). This 
is in good agreement with the value from H09. The distance 
inferred from this distance modulus is 409±23 kpc. 

If we take explicitly into account the dependence of 
Mi.r gbt on global m e tallici ty [M/H], using the calibration 
from Bcllazzini et all (|2004h . we would obtain a distance 
modulus (m-M)o =23.07 and 23.12 for [a/Fe] = 0.0 and = 
+0.4, respectively, fully compatible with the determination 
above . Here we made use of the formula from lSalaris et all 
<ll993h. [M/H] = [Fe/H1+alogl (0.638 xl0 [a;/Fel +0.362) and 
from ba Costa fc Armandrofil l|l990l ). [Fe/H] = -15.16 + 



17.0 x (V - I) 0j _ 3 - 4.9 x (V - I) 2 ,^, where (V - I) 0j _ 3 is 
the mean color of the RGB at absolute I mag = — 3 that for 
our data is (V — I) _ 3 = 1.21. 



5 STELLAR POPULATIONS IN PHOENIX 
AND THEIR SPATIAL VARIATIONS 

5.1 Stellar populations in the Colour-Magnitude 
Diagram 

Figure [5] shows the resulting colour-magnitude diagram 
(CMD) for our FORS mosaic of the Phoenix dwarf galaxy 
(left: Vband; right: /band; lines of Vio ~ 24.8, ho ~ 23.6 
and V$ ~ 25.4, 1$ ~ 24.0 are plotted in the figure). The main 
features of the stellar population of Phoenix are clearly vis- 
ible and indicated with boxes on the CMD: a well defined 
RGB, which should also contain a component of AGB stars; 
a clear detection of an horizontal branch, divided in a red 
(RHB) and blue (BHB) part, extending to V-J~ 0; a well 
populated red clump (RC); an approximately vertical se- 
quence of young main sequence (MS) stars, centered at V- 
/~ —0.2; a vertical sequence of stars emerging from the RC, 
at V-J~0.7, containing bl ue-loop (BL) star s. 

We refer the reader to Holt zman et alJ (|200Ch for a de- 
tailed discussion of the various features. Since we will use 
stars in different evolutionary stages to analyze how the 
stellar population mix varies throughout the galaxy, here 
we focus on how the features present in Phoenix CMD can 
be used as indicators of different age ranges. In evaluating 
the age range dominant for stars in a certain evolutionary 
phase we can use the information on the Phoenix SFH and 
chemical enrichment history derived by H09, i.e. that stars 
older than 6 Gyr have Z between 0.0002 and 0.0004, while 
the stars younger than 2 Gyr have Z between 0.001 and 
0.002. 

From stellar evolutionary models is known that HB will 
contain stars predominantly >10 Gyr old (ancient), while 
the RGB stars will sample the whole stellar population mix, 
with the exception of the stars younger than about 1 Gyr. 
The younger end of the age distribution can be explored 
using the BL and MS stars: the MS stars above the Vio , ho 
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Figure 7. Luminosity function of Phoenix stars (solid histogram) in the magnitude range 18.5 < I< 19.5 (left) and 21.0 < I< 23.2 (right). 
Left: The vertical dotted line indicates the hip, RGB found over about 100 random realizations of the luminosity function convolved with a 
Sobel filter of kernel [-2,0,2] over different binnings (see the text for more details); the dashed line shows one of these random realizations. 
Right: the arrow indicates the position of the RGB bump (see text); the dashed line shows the fit of Eq. (4) to the histogram, while the 
dotted line is the model LF excluding the Gaussian representing the RGB bump. 
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Figure 8. Colour-magnitude diagram for objects classified as stellar in the Phoenix dwarf galaxy (left: V vs V-I ; right: / vs V-I). 
Here we do not include pointing 17, which was observed at a displaced location from the galaxy for checks on the foreground/background 
density. The dashed lines indicate the levels of S/N= 5 and 10 in both bands. The error-bars on the right indicate the typical magnitude 
and color error per magnitude bins of 0.1 mag, and they are derived from the magnitude differences of the stars with double measurements 
(see Fig. [3J . In the left panel we show the selection boxes for stars in different evolutionary phases, as indicated by the labels (see text). 



limit, selected as in [8] are consistent with ages between 0.1- 
0.5 Gyr, while the selected BL stars are sampling slightly 
older stars, mainly 0.5-1 Gyr old. 

The RC contains 1-10 Gyr old stars, in a proportion 
changing with the SFH and metallicity of the stellar popu- 
latio n. Using the predictions on the magnitude of RC stars 
from lGirardi fc Salarisj l|200ll) one can see that for the rather 
narrow and metal-poor metallicity range of stars in Phoenix, 



both the V and I magnitude can act as age indicators. There- 
fore we split the RC stars in magnitude bins of V = 23.2- 
23.4, 23.4-23.6, 23.6-23.8 as indicators of age ranges 2-5 Gyr, 
5-8 Gyr and 8-12 Gyr, respectively. The mean magnitude of 
the RC (derived in Sect. I5.1TT|| indicates that this feature is 
dominated by 5-8 Gyr old stars. 
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5.1.1 Bump on the RGB 

A bump is visible along the RGB and above the RC, both 
in the CMDs in Fig. Hand in the LF (see Fig. [3 right). In 
order to understand whether this bump can be classified as 
an AGB or RGB bump, we derive its mean V and I mag- 
nitudes and compare them to the predictions from stellar 
evolutionary models. 

The I-band LF locally around the RC can be well fit 
with a combination of a polynomial ( for the RGB) and a 
Gaus sian function for the RC stars (jStanek &: Garnavichl 
Il99ct ). To this function we add another Gaussian in order 
to fit the bump present at brighter magnitude than the RC: 

N(I) = a + bxl + cxl+dxe 2 °rc + 

C-Jb,n,p) 2 
g X e 2 "bump ) 

(4) 

where /rc, /bump are the mean magnitude of the RC and 
bump, respectively, and <trc, o"b um p the dispersion. 

We fit Eq. (4) to the LF between 19 < I< 23.5 in bins 
of 0.07 mag in a range of color broadly covering the RC, i.e. 
0.7 - 1.2 (see Fig.0 right). This gives I RC = 22.7 and <ti,rc 
= 0.26 mag for the RC and /bump = 21.90 with (Ti^ump = 
0.044±0.02 mag for the bump. We also fit a similar func- 
tion to the LF in V band over the same colour range and 20 
< V < 24.5; overall, such function gives a good representa- 
tion of the data and yields Vrc = 23.54 with a dispersion 
of <7 V ,rc = 0.22 mag, and Vbump = 22.82 with ov.bump = 
0.13 mag. However, since the bump feature does not seem 
particularly well described by a Gaussian in the Vband, we 
also derive the Vbump from the observed /bump by fitting the 
color distribution of the stars in the range 0.8 < V-I< 1.1 
with a Gaussian. This yields a (V — /)bump = 0.95, result- 
ing in Vbump = 22.85, in good agreement with our other 
determination. 

Using the distance modulus, the extinction and red- 
dening values in Table [1] these would result in absolute 
magnitudes M v ,rc = +0.43, Mi, rc = —0.38, M Vibump = 
—0.26 and Mi >b ump = —1.18. The de-reddened color is 

(V - /)bump,0 = 0.92 

We examine the predicted dependency of the RGB and 
AGB bump magnitude and color as a function of age and 
metallicity, over the range 1 ^ age [Gyr] ^ 12 and 0.0001 
^ Z ^ 0.02, in Fig. I All (top and bottom p anels, respec- 
tively). For this, we use Padua isochrones (|Girardi et al.l 
120001 ; iMarigo et afll20Qg| ). in which the onset and end of 
the various stellar evolutionary phases are conveniently in- 
dicatecQ. We remind the reader that the SFH and chem- 
ical enrichment history derived by H09 for Phoenix show 
that the stars older than 6 Gyr have Z between 0.0002 and 
0.0004, while the stars younger than 2 Gyr old have Z be- 
tween 0.001 and 0.002. Using these isochrones, we find that 
the magnitudes of the detected feature are not compatible 
with those of an AGB bump for the range of metallicities of 
stars in Phoenix and are instead fully compatible the mag- 
nitudes and color of a RGB bump for stars with metallicity 
Z between 0.001 and 0.0001 (see Fig. EH) . 

Since the work of H09 was carried out using Basti 

9 http://stev.oapd.inaf.it/cgi-bin/cmd webpagc 



isochrones (|Pietrinferni et al ] |2004 l 2006) and given that dif- 
ferent sets of isochrones do not always predict similar mag- 
nitudes/colors for stars in the same evolutionary phase (e.g. 
iGallart et alj|2005f ) , we verified if the same conclusion would 
be reached using Basti isochrones. In this case, there is still 
very good agreement between the observed bump magni- 
tudes and colors and those predicted for an RGB bump 
in the metallicity range appropriate for Phoenix stars, but 
there is also a marginal consistency with the values predicted 
for an AGB bump. We then examine the LFs derived for sin- 
gle stellar populations of ages 3-5-8-10 Gyr and Z=0.0001- 
0.001-0.01 using the online tool on the Basti website, and 
we see that the RGB bump is always more populated than 
the AGB bump for the metallicities of Phoenix stars. Since 
in our data we are able to detect only 1 such clump, then 
this is most likely to be the most populated of the two, i.e. 
the RGB bump. Overall, then, the two sets of isochrones 
provide the same conclusion. 

The RGB bump stars sample the intermediate and old 
age populations. For the same age and metallicities the 
AGB bump is expected at similar color, but about 0.5 mag 
brighter in the V-band than the RGB bump; within our 
dataset, however, we do not detect any evident feature at 
that magnitude. 

RGB bumps have been detected in m any LG dwarf 
galax ies, both in dSphs, dlrrs and dTs (see iMonelli et al.l 
I20101 . and references therein); this is the first detection of 
the RGB bump in the Phoenix dwarf galaxy. 

5.2 Spatial distribution of stellar populations 

A view of how the stellar population mix changes across the 
galaxy can already be gleaned by examining the CMD at dif- 
ferent distances from the center. In Figure [9] we plot CMDs 
in bins of projected elliptical radius: R < 2', 2'< R < 4', 4'< 
R < 10'. These CMDs are normalized to have the same num- 
ber of stars, specifically the number of stars above the V5, I5 
limit of the least populated spatial bin (4'< R < 10'). For 
comparison, the CMD at R > 10' in the figure shows where 
the foreground stars and unresolved background galaxies 
should be located. The young MS and BL stars practically 
disappear at R < 4'; the HB (>10 Gyr old) becomes much 
more enhanced going from the inner to the outer parts, and 
the RC becomes less extended in magnitude, indicating a 
decrease in the age range of stars in the stellar population 
mix; the latter is also confirmed by the tendency of the 
mean magnitude of RC stars in the various distance bins 
to become fainter at larger distances from the center, i.e. 
/rc = 22.65, 22.75, 22.86 and V RC = 23.49, 23.60, 23.53 for 
R < 2', 2'< R < 4', 4'< R < 10', respectively. 

Figure [TU] shows the spatial distribution of stars sam- 
pling the different age ranges described in Sect. 15. H i.e. MS, 
BL, RC in bins of V magnitude 23.2-23.4 (RC1), 23.4-23.6 
(RC2), 23.6-23.8 (RC3), and RHB star© We carried out 
a decontamination from foreground stars and background 
galaxies in each of these plots by calculating the predicted 
number of contaminants on the area of the mosaic using the 

10 We plot RHB stars rather than BHB because the latter are in 
smaller number and the selection box may suffer from contami- 
nation from the young MS stars. 
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Figure 9. Colour-magnitude diagram for the objects classified as stellar in the Phoenix dwarf galaxy at different elliptical projected 
radii, R < 2', 2'< R < 4', 4'< R < 10', R > 10' and in the pointing displaced from the region of the mosaic (see labels). To allow for a 
meaningful comparison of the relative strenght of the various features in the CMD, the number of stars in the first 3 panels is normalized, 
so that they contain the same number of stars above the S/N=5 magnitude limit. In the last two panels, whose purpose is only to show 
where the foreground stars and unresolved galaxies are located on the CMD, we do not apply any normalization and show all the stars 
present. The dashed lines indicates the S/N= 5 and 10 limits. The squares w ith error-bars indicate the typical magnitude and color error 
per magnitude bins of 0.1 mag. BASTI isochrones llPietrinferni et al.|[20o3) are overlaid on the CMD (red: 12 Gyr, alpha-enhanced, Z 
= 3 X 10" 4 ; blue: 8 Gyr, solar alpha, Z = 6 X 10 -4 ; the other isochrones have Z = 1 X 10 -3 and ages = 0.1,0.3,0.5,0.7,1.0 Gyr from left 
to righ t). The crosses on the panel R > 10'show the location and amount of Milky Way foreground stars as predicted by the Besancon 
model l lRobin et al.|[2003T) onto the line-of-sight of Phoenix over an area corresponding to the elliptical annuli 10' < R < 13', convolved 
with the observational errors in mag and color; the remaining objects show the location occupied by unresolved background galaxies 
within the observed CMDs. 
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Figure 10. Spatial distribution of stars in different evolutionary stages and age ranges in the Phoenix dwarf galaxy, i.e. MS, BL, RC in 
bins of V magnitude 23.2-23.4 (RC1), 23.4-23.6 (RC2), 23.6-23.8 (RC3), and RHB stars as indicated by the labels and boxes in Fig. GO 
The various populations have been extracted from the S/N>10 sample. A decontamination has been carried out, using the contaminant 
density derived as a weighted mean of the outer points of the surface density profile (except for the MS stars for which we used the 
contaminant density derived from the objects at R > 12.6'). The ellipse shows the measurement of the nominal tidal radius derived in 
this work. A red cross is placed at the coordinates of Phoenix optical center to guide the eye. North is up and East is on the left. 
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contaminant density derived from method a) applied to all 
these individual stellar populations, except for the MS stars 
to which we applied method b) because of their asymmetric 
distribution clearly different from the rest of the population 
(see the figure); in practise, the number of contaminants 
was randomly removed from the objects falling in a specific 
selection box, taking care of the fact that the density of 
contaminants is mostly expected to be uniform over such a 
small area. 

It is directly visible from Figure [10] that the younger 
the stellar population, the more centrally concentrated is 
its spatial distribution. This is particularly evident for the 
youngest stars, the MS and the BL, which also show a less 
extended distribution with respect to the rest. This would 
be consistent with the star formation region shrinking with 
time, with recent star formation having occurred in the in- 
nermost regions where the gas density would have still been 
high enough to sustain it (see also H09). 

Such spatial variations of the stellar population mix 
had also been observed in previous studies usin g shallower 
and / or less spatially extended data (e.g. MD99, Held ct al. 
1 19991 . H09). The depth and spatial extent of our data allows 
us to perform a more accurate analysis than previous stud- 
ies and quantify the variations in the spatial distribution of 
stars by deriving a surface number density profile for each 
population covering a different evolutionary stage (except 
for the flattened, asymmetric MS stars). The methodology 
and the functional forms fitted to the observed profiles are 
as described in Sect. 13.21 The best-fitting parameters (see 
Table [7]) show that the variations among the spatial dis- 
tribution of the various stellar populations appear smooth, 
slowly and steadily going from more concentrated distri- 
butions from the youngest populations to more extended 
for the oldest. We checked that using the contaminant den- 
sity derived from method b) to all stellar populations would 
bring negligible changes to the resulting spatial distribu- 
tions, surface brightness profiles and resulting best-fit pa- 
rameters (as in the case of the overall stellar population, 
the reduced % 2 values of the best-fitting profiles are larger 
when using method b) than method a) ). Also, the deter- 
mination of the outer points over which we calculate the 
density of contaminants was adjusted by identifying where 
the surface number profile was becoming approximately flat; 
we checked that the choice of this region did not influence 
the results. 

As noted in previous works (e.g. MD99), the MS stars, 
those approximately younger than 0.7 Gyr, have a flattened 
distribution, almost perpendicular to the other stellar com- 
ponents. There appears to be two main clumps of young MS 
stars. Some clumping is also present in the distribution of the 
other young stars - the BL and the RC1 that represent the 
stars with ages mostly between 2-5 Gyr -, which both show a 
(rather spherical) concentration of stars displaced from the 
center, towards the South-East at rectangular coordinates 
(+0.5,-0.5); this is the direction of the coordinates we find 
for the position of the center in the inner 2' in Fig. [5] Such 
displaced concentration disappears when moving to older 
ages, where the distribution becomes more regular. It can 
be noted that the populations on average older than 2 Gyr 
have a regular, spheroidal appearance and do not display a 
disk-like distribution as the young MS stars. Therefore the 
overall configuration of Phoenix does not resemble a disk- 



halo like larger spiral galaxies such as the MW: in the MW 
the disk contains stars of all ages (from young to ancient) 
while the halo is formed solely of ancient stars; this is differ- 
ent from the situation in Phoenix, where the disk-like feature 
is composed only of young stars, while the spheroidal part of 
the galaxy contains stars from intermediate to ancient ages. 

We can speculate that the flattened, asymmetric dis- 
tribution of the young MS stars still retains the imprint of 
where the recent star formation took place, and therefore of 
where the gas was recently located. This appears also to be 
still imprinted in the distribution of slightly older stars, such 
as the BL and the RC1, in the form of the clump displaced 
from the center. However, being slightly older than the MS, 
these BL and RC1 stars would have had more time to be- 
come sensitive to the overall potential of the galaxy, and for 
their spatial distribution to become smoother. 



6 DISCUSSION 

• Spatial variations of the stellar population mix: The spatial 
distribution of the stars changes with radius and becomes 
less and less centrally concentrated the older the stellar pop- 
ulation. This is reflected in the values of the best-fitting pa- 
rameters of the surface number profiles of the different stel- 
lar populations that steadily go from a more concentrated 
to less concentrated distribution with increasing ages. 

Similar variations in the stellar population mix appear 
to be a common characteristic of LG dwarf galaxies (see 
e.g.lHarbeck et al.ll200ll; iTolstov et al.ll2004l;lBattaglia et all 



l2006l ; lBernard et al.ll2008l ; lMonelli et al.ll2012l ~ with the star 



formation proceeding longer in the central regions. There 
are a number of mechanisms that could cause this: a nat- 
ural evolution of the gas, for example p rogressively sinking 
in the center in abse nce of rotation (e.g. IStinson et al.ll2009l ; 
ISchroven et al.|[201 lh or reaching higher densities in the cen- 
ter and therefore continuing star formation for longer. An- 
other explanation could be ejection of the ISM because of su- 
pernovae explosions being more efficient in the outer parts; 
removal of the ga s from ram pressure/tidal stripping (e.g. 



iMaver et al1l200fj 



Stinson et al.l ((2009) show that age gradients of this 



type can arise in model isolated galaxies of low mass (with 
virial velocities <20 km s _1 ) without any need for exter- 
nal factors, because these systems do not form stable star 
forming discs, as the gas is mostly supported by pressure 
rather than angular momentum. Star formation rate is ini- 
tially high in the center, because of the larger densities of the 
gas, and while the gas is consumed, the disc contracts until 
pressure support is re-established. Further shrinking may be 
due to a lower turbulence velocity due to the decline in SFR 
in the center, which mantains high gas densities in the cen- 
tre but not in the outer regions, resulting in star formation 
that is systematically less radially extended as time goes on. 

It appears that whatever the mechanism is that creates 
evident variations in the stellar population mix, it can act on 
very different timescales from galaxy to galaxy. The Sculp- 
tor dSph, w hich predominantly c onsists of stars older than 
10 Gyr (e.g. Ide Boer et aLll2012T ). clear ly displays this fea- 
ture already in its ancient population (|Tolstov et al.l [20041 ; 
Ide Boer et al.ll2~012l ). In the Fornax dSph instead there is no 
detected spatial variation in the mix of the ancient popula- 
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tion, while these changes are clear when comparing ancient , 
intermediate age and young stars (e.g. IStetson et al.lll99sl; 
iBattadia et alJ l2006). Phoenix in this respect is similar to 
Fornax. Such spatial variations between populations of dif- 
ferent ages could be expected if they are due to environ- 
mental effects: more isolated objects and satellites with less 
internal/eccentric orbits around the host galaxy would be 
less affected and star formation could proceed longer in the 
outer parts; this would be consistent with what is presently 
known about the orbits of Sculptor and Fornax, bar the large 
error-bars in the proper motions values. At present there are 
no determinations of the proper motion of Phoenix; assum- 
ing a radial orbit for this object, its systemic velocity (e.g. 
Ilrwin Sz Tolstoy||2002l ) referred to the Galactic Rest Frame 
(vsys.GSR ~ —100 km s _1 ) would yield an r x v 2 about 1.5 
times lower than for Leo I, which is considered as a satellite 
of the Milky Way: in principle Phoenix could therefore still 
be a Milky Way satellite falling back. 

In order to disentangle environmental effects from inter- 
nal ones, it would be instructive to quantify the differences 
in the distribution of stellar populations of various ages in 
a sample of LG dwarf galaxies. The isolated objects would 
give insights on the extent to which internal mechanisms 
can be responsible for the variations in the stellar popula- 
tion mix. It would be also instructive to compare systems 
with similar SFHs: for example Sculptor, Tucana and Ce- 
tus are all dSphs which produced most of their stars very 
early on, but while the former is a MW satellite, the other 
two are found far from the MW and M31. There are indica- 
tions that also among the ancient stars in Tucana and Cetus, 
the older/metal-poor ones have a more extended spatial dis- 
tribution than the younger /metal-rich (|Harbeck et al.l[200ll : 
iBernard et~afl 120081 : [Monelli et all 120121) ; it remains to be 
quantified how these differences in the distribution of the an- 
cient stars compare between isolated and non-isolated dwarf 
galaxies. 

• Spatial distribution of young stars: The stars younger 
than 0.7 Gyr display a clearly different spatial distribution 
than the rest of the population: not only they virtually dis- 
appear at distances larger than 4' from the center, but they 
are almost perpendicular to the bulk of the other stars, and 
the distribution is asymmetric and more flattened. This was 
already noted by VDK91 and MD99. By separating the sam- 
ple in age bins, we also see that some asymmetries in the 
form of clumps are visible in the spatial distribution of BL 
stars and the younger RC stars, while stars older than about 
5 Gyr display regular spatial distributions. 

We can speculate that the flattened, asymmetric dis- 
tribution of the young MS stars still retains the imprint 
of where recent star formation took place, and therefore of 
where the gas was recently located. This appears also to be 
still imprinted in the distribution of the BL and young RC 
stars; however, since these stars are older than the young MS 
stars, they have had more time to diffuse to larger distance 
and to become more sensitive to the overall potential of the 
galaxy, resulting in a smoother distribution. Following this 
reasoning, it is likely that in a few 100s Myr also the young 
MS stars will diffuse to larger distances and acquire a more 
regular morphology, similar to the one of the main body of 
the galaxy. At that stage, Phoenix would probably look like 
a typical dSph, if it also were to lose its HI gas or exhaust 
it in further star formation. 



In the models of Stinson et al. (2009) for low mass ob- 
jects, although the majority of stars form systematically 
more inwards the more time passes, some amount of ra- 
dial migration is also present, which goes in the direction 
of pushing stars to larger radii with respect to their birth 
site, qualitatively consistent with what we see here for the 
youngest stars. 

We point out that a feature such as the one seen here 
- an inner flattened component containing young stars and 
tilted with respect to the main body of the galaxy - is also 
present in ot her LG dwarfs su ch as the Leo A dlrr (e.g. 
see Fig. 1 in ICole et al l I2007T ) and the closer by Fornax 
dSph, a MW satellite. With respect to the majority of MW 
dSphs, Fornax has had a much more extended SFH and is 
the only one to show star s possibly as young as 50 Myr 
(I Coleman fc de Jond 120081 ). Also in this galaxy the young 
MS stars show a flattened distribution, tilted with respect 
to the main stellar population, while the slightly older BL 
stars have a less asymmetric and flattened distribution and 
appear to have diffused to larger distances. In this case, the 
tilting of the y oung stars is less e nhanced, being approxi- 
mately 40° (e.g. IStet son et all 1 19981 ). 

It is unclear what would cause the tilting of the youngest 
stars with respect to the bulk of the stellar population. 
Such a behaviour has not been reported for more mas- 
sive, clearly rotating syste ms, such as for example the WLM 
dlrr (jLeaman et al.ll2012l ). where the young stars would be 
aligned and concentric with the distribution of older stars 
and the gas. This may indicate an influence of the total 
mass of the galaxy in determining its intrinsic angular mo- 
mentum, and consequently the resulting distribution of its 
stellar populations: we speculate that, in less massive sys- 
tems, the low amount of angular momentum would permit 
a disordered configuration of the gas, and therefore of the 
stars young enough not to have become sensitive to the 
overall potential of the galaxy. In apparent contradiction 
to that, in NGC 6822, a dlrr with a luminosity similar to 
WLM, the spatial distribution of the young stars and the old 
RGB stars are almost perpendicular to each other, with th e 
young stars following the HI disk fe.g. iDemers et ai1l2006l ): 
in this case however the picture is comp licated by the com- 
plex disk structure (|Cannon et al.l |2012| ) and the classifica- 
tion of NGC 6822 as a polar ring galaxy (e.g. IDemers et al.l 
l2006h . 

• Disk-halo appearance: Previous works (e.g. MD99) 
showed that the central parts of Phoenix display an inner 
component that is tilted with respect to the main body of 
the galaxy and that morphologically resembles a disk - while 
the main body has a spheroidal distribution. The question 
was raised whether Phoenix may have a disk-halo structure 
similar to that of large spirals galaxies, such as the MW. 

We can comment on this issue simply making some con- 
sideration on the age of stars found in the "disk" and "halo" 
of Phoenix. 

In this work we showed that in Phoenix only young stars 
are responsible for the elongated central structure that mor- 
phologically resembles a disk, while such structure is absent 
for stars older than 2 Gyr, which display a rather regular 
spheroidal morphology. Therefore only stars younger than 
1 Gyr are found in the "disk", while both intermediate age 
and ancient stars are found in the "halo". This is clearly 
different to what is observed for the MW, where the disk is 
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known to contain stars of all ages while the halo consists of 
ancient stars; therefore the structure we see in Phoenix is 
not analogous to what found for the MW and M31. 

This result is also difficult to reconcile with the prop- 
erties of stellar haloes formed in a ACDM context. Stellar 
haloes formed from disrupted satellite galaxies are a natural 
consequence of the hierarchical formation of galaxies in a 
ACDM framework and, for large spirals, the average age of 
stars in the stellar haloes is expected to be around 11 Gyr 
l|Cooper et al.ll2010h . with essentially no stars younger than 
5 Gyr, from mainly a few satellites accreted between redshift 
1 and 3. Since in this context dwarf galaxies are expected 
to be among the first galactic systems to form, their stellar 
haloes should arguably have formed at earlier times than 
for large spirals, and from smaller progenitors, and should 
therefore contain only ancient stars. 



7 SUMMARY AND CONCLUSIONS 

We presented results from wide-area photometry in Fand 
/ band from VLT/FORS for Phoenix, one of the few tran- 
sition type dwarf galaxies of the LG. The data consist of a 
mosaic of images covering 26' x 26'around the optical cen- 
ter of the galaxy and reaching below the horizontal branch 
of the system: this is the only data-set for this galaxy that 
combines relatively deep photometry with such a large spa- 
tial coverage. For comparison with previous studies, we can 
trace Phoenix stellar population out to projected radii of 1.5 
kpc from its center, while the study of H09 reached out to 
0.5kpc. 

One output of the analysis is the re-determination of the 
overall structure of the system and especially of a determina- 
tion of its extent, that was very uncertain in the literature. 
We derived a nominal tidal radius of r± = 10.56' ± 0.15', 
about 60% the value determined by MD99. The best-fitting 
profile to the overall stellar population of Phoenix is a Sersic 
profile of Sersic radius Rs = 1.82'±0.06' and m =0.83±0.03. 

The distance of Phoenix has been derived indepen- 
dently from our dataset using the RGB tip method, and it 
is in agreement with previous distance determinations from 
H09 and M99. Examination of the CMD and LF revealed 
the presence of a bump above the red clump, consistent with 
being a RGB bump. This is the first detection of this feature 
in Phoenix. 

The depth of the photometry allows to study stars in 
different evolutionary phases, with ages ranging from about 
0.1 Gyr to the oldest ages as they can be identified from fea- 
tures in the CMD. This, combined to the large area covered, 
enabled us to explore in a very direct way how the stellar 
population mix varies across the face of the galaxy. 

The spatial distribution of the stars changes with radius 
and becomes less and less centrally concentrated the older 
the stellar population. This is reflected in the values of the 
best-fitting parameters of the surface number profiles of the 
different stellar populations that we provide. 

As also found in previous studies, Phoenix displays in 
the inner regions a flattened structure, almost perpendicular 
to the main body of the galaxy. This feature was suggestive 
of a disk-halo structure and was raising the question whether 
also Phoenix could consist of a a disk-halo system such as the 
one found in the Milky Way. We showed that the flattened, 



tilted structure appears to be present mainly among stars 
younger than 1 Gyr, and absent for the stars > 5 Gyr old, 
which on the other hand show a regular distribution also 
in the center of the galaxy. This argues against a disk-halo 
structure of the type found in large spirals such as the Milky 
Way. 
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APPENDIX A: RGB AND AGB BUMP 

This paper has been typeset from a TgX/ DTfjX file prepared 
by the author. 
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Figure Al. For this figure we used solar-scaled Padua isochrones (Marigo et al. 2008) spanning an age range from 2 to 13 Gyr, Z=0.0001- 
0.02. Top: The symbols connected by vertical lines show the luminosity interval of RGB stars crossing the discontinuity in chemical profile 
left by the first dredge-up event. This interval corresponds to the bump in the luminosity function along the RGB. Bottom: as before but 
this interval corresponds to the clump of early-AGB stars in colour-magnitude diagrams, therefore the AGB bump. From left to right 
we show how the V mag, I mag and V-Ivaxy as a function of age in Gyr for isochrones of different metallicities. 



